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Carbon  disulfide  is  the  most  popular  material  for  applications  of  nonlinear  optical  (NLO)  liquids,  and  is 
frequently  used  as  a  reference  standard  for  NLO  measurements.  Although  it  has  been  the  subject  of  many 
investigations,  determination  of  the  third-order  optical  nonlinearity  of  CS2  has  been  incomplete.  This  is  in 
part  because  of  several  strong  mechanisms  for  nonlinear  refraction  (NLR),  leading  to  a  complex  pulse  width 
dependence.  We  expand  upon  the  recently  developed  beam  deflection  technique,  which  we  apply,  along 
with  degenerate  four-wave  mixing  and  Z-scan,  to  quantitatively  characterize  (in  detail)  the  NLO  response  of 
CS2,  over  a  broad  temporal  range,  spanning  6  orders  of  magnitude  (~32  fs  to  17  ns).  The  third-order 
response  function,  consisting  of  both  nearly  instantaneous  bound-electronic  and  noninstantaneous 
nuclear  contributions,  along  with  the  polarization  and  wavelength  dependence  from  390  to  1550  nm, 
is  extracted  from  these  measurements.  This  paper  provides  a  self-consistent,  quantitative  picture  of  the 
third-order  NLO  response  of  liquid  CS2,  establishing  it  as  an  accurate  reference  material  over  this  broad 
temporal  and  spectral  range.  These  results  allow  prediction  of  the  outcome  of  any  NLR  experiment 
on  CS2.  ©  2014  Optical  Society  of  America 

OC/S  codes:  (190.4400)  Nonlinear  optics,  materials;  (300.6420)  Spectroscopy,  nonlinear;  (190.3270)  Kerr  effect;  (320.7110)  Ultrafast 
nonlinear  optics. 
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1.  introduction 

Carbon  disulfide  (CS2)  is  a  widely  used  nonlinear  optical 
(NLO)  liquid  owing  to  its  large  third-order  nonlinear  refrac¬ 
tion  (NLR).  It  has  been  subject  to  many  experimental  studies 
utilizing  time-resolved  techniques  such  as  optical  Kerr  effect 
(OKE)  [J_-8],  degenerate  four-wave  mixing  (DFWM) 
[9,10],  and  nonlinear  interferometry  [1 1],  as  well  as  frequency 
domain  light  scattering  [12—141,  third-harmonic  generation 
[15],  and  Z-scan  [16-211.  It  is  found  in  a  wide  array  of 


NLO  applications,  including  liquid-core  optical  fibers  for  non¬ 
linear  photonics  applications  [22,23],  soliton  propagation 
[24,25],  supercontinuum  generation  [26],  slow  light  [27],  ul¬ 
trafast  time- resolved  imaging  [28],  and  all-optical  switching 
[29] .  In  addition,  it  is  commonly  used  as  a  reference  material 
for  nonlinear  spectroscopy  [15,30,311 .  However,  the  NLO  re¬ 
sponse  of  CS2  varies  by  more  than  1  order  of  magnitude,  over 
the  time  range  of  ~30  fs  to  ~30  ps  [17].  This  is  because  CS2 
exhibits  a  strong  noninstantaneous  third-order  response, 
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arising  from  motions  of  the  molecules.  Two-photon  absorp¬ 
tion  (2PA)  also  begins  to  play  an  important  role  in  the  non¬ 
linear  response  for  wavelengths  shorter  than  500  nm  [321,  and 
via  Kramers-Kronig  relations  the  bound-electronic  compo¬ 
nent  of  the  NLR  displays  dispersion  in  this  region  [331-  With¬ 
out  broad  and  accurate  characterization  of  the  pulse  width  and 
wavelength  dependence  of  its  nonlinearity,  the  use  of  CS2  as  a 
calibration  material  is  problematic. 

Many  of  the  previous  time-domain  measurements  of  CS2 
have  focused  solely  on  relative  magnitudes  (and  dynamics)  of 
the  nuclear  response,  and  provide  no  information  on  the  ab¬ 
solute  magnitude  of  the  nonlinearity  [J_, 5,6,9, J_0,34] .  Other 
time  domain  investigations  have  used  a  relative  measurement 
technique,  and  compared  their  results  to  literature  values  of  the 
nonlinearity  of  a  reference  material  [JJJ.  Reports  of  quantita¬ 
tive  measurements  of  NLR,  using  absolutely  calibrated  tech¬ 
niques  (such  as  Z-scan),  only  pertain  to  specific  pulse 
widths  and  wavelengths,  and  are  thus  limited  in  applicability. 

In  this  paper,  we  present  a  thorough  experimental  investi¬ 
gation  of  the  third-order  NLO  response  of  CS2.  By  applying 
our  newly  developed  beam  deflection  (BD)  technique  [351,  we 
measure  the  third-order  temporal  response  function,  including 
the  absolute  magnitude  and  symmetry  of  each  component, 
without  a  need  to  scale  to  a  reference  material.  The  BD 
method  is  particularly  useful,  for  measuring  the  absolute  mag¬ 
nitudes  of  the  various  nonlinear  mechanisms,  by  displaying  the 
temporal  dynamics  and  polarization  dependence.  We  then 
apply  the  measured  response  function  to  predict  the  results 
of  additional  independent  NLO  experiments:  DFWM  and 
Z-scans,  with  pulses  from  32  fs  to  17  ns,  including  linear 
and  circular  polarizations.  Additionally,  we  determine  the 
dispersion  of  NLR,  including  both  bound-electronic  and 
nuclear  contributions,  and  the  spectrum  of  2PA  from  390 
to  1550  nm.  These  data  give  a  self-consistent,  quantitative 
picture  of  the  third-order  nonlinear  dynamics  of  liquid  CS2, 
making  possible  a  prediction  of  the  outcome  of  any  measure¬ 
ment  of  2PA  or  NLR. 


the  molecular  ensemble.  Each  contribution  to  the  third-order 
response  causes  an  irradiance,  and  time-dependent  change  in 
the  refractive  index,  which  for  a  single  beam  is  given  by 


A n(t)  =  n2ie/I(t)  +  /  R(t  - 


(1) 


where  I(t)  is  the  irradiance,  R(t)  is  the  noninstantaneous  com¬ 
ponent  of  the  third-order  response,  and  «2,r/  >s  the  nonlinear 
refractive  index  that  originates  from  the  bound  electrons.  All 
experiments  are  outside  the  linear  absorption  regime,  and  thus 


_3Re(/(3)) 
ni,d  4«q£0C 


(2) 


where  »0  is  the  linear  refractive  index,  e0  the  vacuum  permit¬ 
tivity,  the  bound-electronic  third-order  susceptibility,  and 
c  is  the  speed  of  light  in  vacuum  [15].  The  first  term  of  Eq.  (1) 
describes  the  bound-electronic  response  and  follows  the  tem¬ 
poral  irradiance  profile  of  the  pulse,  while  the  second  term 
describes  the  change  in  refractive  index  due  to  the  noninstan¬ 
taneous  nuclear  motions.  While  there  have  been  theoretical 
investigations  of  the  effect  of  correlations  of  various  nuclear 
motions  on  the  third-order  optical  response  [36,371,  multiple 
investigators  have  found  good  agreement  with  experiments  by 
treating  these  mechanisms  as  linearly  independent  [1,5,1 1], 
We  take  the  same  approach  here,  and  write  the  nuclear  con¬ 
tribution  as  a  sum  of  three  responses: 

R(t)  =  S>2n2,mrm(t),  (3) 

m 


where  «2  m  is  the  magnitude  of  the  wzth  mechanism  and  the 
temporal  response  function  rm{t)  is  normalized  such  that 


M&t  =  1. 


(4) 


2.  ORIGINS  OF  NONLINEAR  REFRACTION  IN  CS2 

Several  mechanisms  give  rise  to  a  change  in  the  refractive  index 
of  CS2,  each  with  its  own  characteristic  response  time  and 
polarization  dependence.  Within  the  Born— Oppenheimer 
approximation,  these  mechanisms  may  be  divided  into  two 
categories:  those  of  the  bound  electrons,  and  those  of  the 
nuclei.  The  bound-electronic  contribution  to  the  third-order 
response  originates  from  the  purely  electronic  second  hyperpo¬ 
larizability  of  the  individual  molecules.  The  response  time  is 
limited  by  the  time-energy  uncertainty  relation,  to  typically 
<  1  fs,  which  is  essentially  instantaneous  compared  to  the  pulse 
duration,  and  is  related  to  the  transition  dipole  coupling 
between  discrete  electronic  states  of  the  molecule.  The  char¬ 
acteristic  response  times  of  the  nuclear  mechanisms  to  the 
index  change  are  approximately  hundreds  of  fs,  owing  to 
the  large  mass  of  the  nuclei.  These  nuclear  motions  give  rise 
to  a  third-order  optical  response,  by  altering  the  bulk  linear 
susceptibility  either  through  a  change  in  the  molecular 
polarizability,  or  by  changes  in  the  average  orientation  of 


The  dominant  nuclear  mechanism  in  CS2  is  molecular  re¬ 
orientation  [38] .  Due  to  its  linear  shape  it  exhibits  anisotropic 
polarizability  and  has  no  permanent  dipole  moment.  A  strong 
electric  field  induces  a  dipole  leading  to  an  applied  torque, 
causing  the  molecular  axis  to  align  with  the  incident  field. 
The  time  scale  of  reorientation  depends  on  the  moment  of 
inertia  of  the  molecules,  which  are  angularly  accelerated  by 
the  electric  field.  After  a  sufficiently  short  pulse,  the  molecules 
will  continue  to  orient  toward  the  field  and  then  thermally 
randomize  at  a  rate  that  depends  on  the  viscosity  of  the 
liquid.  Since  the  diffusive  nature  of  the  relaxation  follows 
the  Debye-Stokes-Einstein  relation  [39],  this  mechanism  is 
referred  to  as  diffusive  reorientation.  The  molecules  may  be 
modeled  as  a  driven  damped  harmonic  oscillator  [JJ,  where 
in  the  overdamped  case  the  temporal  response  function  is 
well-described  by 

rd(t)  =  Cd[\  -  e  r^e  Tf-d&(t)  , 


(5) 
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where  the  subscript  d  indicates  the  diffusive  reorientation 
mechanism,  r rj  and  ^  are  the  rise  and  fall  times,  respec¬ 
tively,  C d  =  (rrij  +  t fj)  / 'tjr  d  is  a  normalization  factor, 
and  0(r)  is  the  Heaviside  function  which  ensures  causality. 

In  liquids,  the  interactions  of  neighboring  molecules 
must  be  considered  when  looking  at  molecular  motion. 
One  example  is  libration,  which  is  an  oscillatory  rocking 
motion  resulting  from  a  molecule  aligning  toward  the  electric 
field  but  hindered  by  the  environment.  The  influence  of  this 
mechanism  on  the  refractive  index  is  due  to  the  coherent 
oscillation  of  many  molecules.  Hence,  a  loss  of  coherence 
between  many  oscillators,  caused  by  collisional  dephasing,  is 
the  dominant  decay  mechanism  of  this  induced  index  change 
[J_] .  The  quantum  harmonic  oscillator  model  of  the  (under¬ 
damped)  librational  response  function  is  given  by  [4] 

rl{i)  =  Cle  £©(*)  (6) 

Jo  CO 

where  the  subscript  /  indicates  the  librational  mechanism  and 
g(co)  is  the  distribution  function  of  libration  frequencies.  Since 
different  molecules  experience  different  environments  within 
the  liquid,  this  librational  motion  is  inhomogeneously  broad¬ 
ened  [40].  We  use  the  “antisymmetrized”  Gaussian  distribu¬ 
tion  function  [14] 

_  ( co-cop )2  _  (to+co o)2 

g(o>)  =  e  z.2  -  e  i"2  ,  (7) 

where  tu0  is  the  central  frequency  and  a  is  the  spectral  width. 

Finally,  there  are  intermolecular  collision-induced  varia¬ 
tions  in  the  molecular  polarizability.  One  cause  is  dipole- 
induced  dipole  interactions,  where  in  addition  to  the  applied 
field,  each  molecule  also  experiences  the  reradiated  electric 
field  from  the  induced  dipole  in  the  neighboring  molecules, 
which  distorts  its  polarizability  [41] .  Electronic  overlap  effects 
become  important  for  very  small  molecular  separations 
[37,41] .  Time  domain  measurements  have  found  that  the  col¬ 
lision-induced  mechanism  follows  an  exponential  rise  and 
decay  function  [J_, JT  ,42] .  The  (overdamped)  collision-induced 
contribution  follows  a  response  function  of  the  form 

r,(*)  =  C,(  1 -,-£)/£©(*),  (8) 

where  the  subscript  c  indicates  the  collision-induced  mecha¬ 
nism.  Similar  responses  have  been  observed  in  atomic  liquids, 
including  Ar  and  Xe  [42],  which  have  isotropic  polarizabilities. 

The  four  nonlinear  mechanisms  have  two  different  polari¬ 
zation  dependencies.  Given  the  isotropic  nature  of  the  liquid  at 
room  temperature,  the  bound-electronic  and  collision-induced 
mechanisms  follow  isotropic  symmetry.  However,  diffusive  re¬ 
orientation  and  libration  have  a  quite  different  polarization 
dependence.  As  the  CS2  molecules  align  toward  the  electric 
field,  anisotropy  is  generated  within  the  liquid.  The  average 
polarizability,  and  thus  refractive  index,  is  increased  in  the 
direction  of  alignment  and  decreased  in  the  perpendicular 
directions.  We  refer  to  this  type  of  polarization  dependence 
as  reorientational.  The  isotropic  and  reorientational 


mechanisms  may  be  easily  isolated  by  conducting  measure¬ 
ments  with  the  appropriate  polarizations.  This  gives  an 
additional  dimension  to  aid  in  the  characterization  of  the 
third-order  response  of  CS2. 

3.  BEAM  DEFLECTION 

To  fully  characterize  the  temporal  dynamics  of  CS2,  we  make 
use  of  our  recently  developed  BD  technique  [35].  As  shown  in 
Fig.  J_,  a  strong  excitation  beam  induces  a  change  in  the  refrac¬ 
tive  index  following  its  Gaussian  spatial  profile.  The  resulting 
index  gradient  is  sensed  by  a  weaker,  and  spatially  smaller, 
probe  beam  aligned  off  center  of  the  excitation  beam  where 
the  irradiance  gradient  is  maximized.  Typically,  the  probe 
beam  has  a  size  ~5x  smaller  than  that  of  the  excitation  beam. 
The  prism-like  index  gradient  induced  by  the  excitation  causes 
the  probe  to  be  deflected,  as  measured  by  a  segmented  quad¬ 
cell  detector  in  the  far-field.  The  BD  technique  has  distinct 
advantages  over  other  nonlinear  excite-probe  techniques.  It 
produces  a  signal  that  is  linearly  proportional  to  the  induced 
refractive  index  change,  directly  giving  its  sign,  and  thus  does 
not  require  optical  heterodyne  detection  as  OKE  or  DFWM. 
In  addition,  various  tensor  elements  may  be  measured  directly 
by  altering  the  polarizations  of  the  excitation  and  probe  beams, 
whereas  OKE  measures  the  induced  birefringence.  BD  is  much 
easier  to  implement  than  two-color  Z-scan  [43]  or  interfero¬ 
metric  techniques  [1 1], 

To  eliminate  experimental  complications,  such  as  coherent- 
coupling  effects  [5,44,45]  and  a  large  background  signal  from 
the  excitation,  we  use  different  wavelengths  for  the  excitation 
and  probe.  However,  the  relatively  high  dispersion  of  CS2  [46] 
causes  the  pulses  to  propagate  with  different  group  velocities 
and  pass  through  one  another  in  time,  resulting  in  a  reduced 
signal  [47].  Following  Negres  et  al.  [47],  we  develop  a  simple 
method  of  analysis  for  the  BD  measurements  in  the  presence  of 
group  velocity  mismatch  (GVM)  and  a  noninstantaneous 
nonlinear  material  response.  This  may  be  applied  to  any 
material  in  the  limit  where  the  pulse  widths  do  not  change 
upon  propagation  through  the  sample,  i.e.,  negligible  group 
velocity  dispersion  (GVD).  In  our  experiments,  pulses  are  es¬ 
timated  to  broaden  by  less  than  9%  [46] .  This  analysis  allows 
us  to  experimentally  determine  the  third-order  nonlinear 
response  of  CS2,  including  absolute  magnitudes  of  each 


Fig.  1.  (a)  Schematic  diagram  of  BD  experiment  [35];  (b)  irradiance 

distribution  of  excitation  (red),  and  probe  (green)  at  sample  plane; 
(c)  zoom  in  of  (a)  at  sample  showing  a  profile  of  A n\  (d)  probe  beam 
on  segmented  detector  without  and  with  deflection. 
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mechanism,  their  characteristic  response  times,  and  polariza¬ 
tion  dependence. 

Equation  (9)  gives  the  nonlinear  wave  equation  for  a  weak 
probe  centered  at  frequency  go  in  the  limit  of  zero  depletion  of 
the  excitation  pulse  at  coe  [48].  The  sample  is  assumed  to  be 
thin,  meaning  the  thickness  is  less  than  the  Rayleigh  range  of 
both  beams,  and  the  peak  induced  phase  changes  are  small 
[161.  This  makes  the  excitation  irradiance  Ie  independent 
of  propagation  distance  z  and  the  deflection  of  the  probe 
within  the  sample  negligible.  The  probe  electric  field  is 
E(r ,z,t)  =  1  /28,{r,  z,  fee,  where  £(r,z,t)  is 

the  field  amplitude,  r  is  the  radial  coordinate,  k  =  aon0/c 
is  the  wave  number  and  e  is  the  electric  field  unit  vector. 

t(jfz  +  l£)£  +  ko£(l nXelIe(t)  +  /“  R(t  -  f')/,(Odf') 

=  0,  (9) 

where  k0  =  co/c  and  v  is  the  probe  group  velocity  (see 
Supplement  1  for  derivation).  The  second  term  of  Eq.  (9)  rep¬ 
resents  the  nonlinear  interaction  of  the  probe  pulse  with  the 
excitation.  It  is  written  as  the  sum  of  the  bound-electronic  con¬ 
tribution  (plus  the  convolution  of  the  excitation  pulse),  with 
the  noninstantaneous  nuclear  component  of  the  third-order 
response  function.  The  factor  within  the  parentheses  is 
A n(t),  experienced  by  the  weak  probe,  induced  by  the  exci¬ 
tation.  Note  the  instantaneous  component,  2n2e{I  e(t) ,  is 
greater  by  a  factor  of  2  than  that  in  Eq.  (J_).  The  excitation 
and  probe  fields  interfere,  and  generate  a  refractive  index  gra¬ 
ting,  off  which  the  excitation  diffracts  into  the  probe  direction 
(and  vice  versa),  with  the  appropriate  phase  and  amplitude  to 
double  the  phase  shift,  so  called  weak-wave  retardation 
[44,49].  This  is  true  even  in  the  nondegenerate  case,  since 
the  instantaneous  bound-electronic  response  is  able  to  follow 
the  rapidly  moving  grating.  The  noninstantaneous  nuclear 
component,  however,  is  too  slow  to  follow  the  grating,  and 
no  additional  phase  shift  appears  [44] .  In  the  degenerate  case, 
this  grating  is  stationary,  or  very  slowly  moving  for  chirped 
pulses,  and  the  noninstantaneous  response  does  contribute  ad¬ 
ditional  terms  to  the  nonlinear  interaction  [45].  Equation  (9) 
assumes  the  nondegenerate  case  Goe  ^  go. 

Following  [47],  the  solution  for  Gaussian  pulses  with  the 
addition  of  the  noninstantaneous  third-order  response,  and 
the  transverse  spatial  profiles  of  the  excitation  and  probe,  is 
given  by 


cv„  t a  c  (  (T  +Td-p)2  ,k0L 

£(r,  T)  =  c0(r)  exp  I - ^ - \-i—IeS>(r) 

■  1 2«2,<./[erf (T)  -  erf (7  -  p)\ 

+  JJ  J°°r(t2-  7’1K7idr1dr2|)J  (io) 

where  T  =  (t  -  z/ve)  /re  is  the  dimensionless  normalized 
time  (that  moves  with  the  excitation  group  velocity  ve), 
re  is  the  excitation  pulse  width  (HWl/eM),  Td  is  the 


dimensionless  delay  between  excitation  and  probe,  T  = 
r / re  is  the  normalized  probe  pulse  width,  and  L  is  the  sample 
thickness.  The  GVM  parameter  is 


Te\V  Vj 


ng,e)> 


(ID 


where  ng  and  nge  are  group  indices  of  the  probe  and  excitation, 
respectively.  Equation  (10)  describes  the  propagation  of  a  weak 
probe  pulse  accounting  for  both  instantaneous  and  noninstan¬ 
taneous  nonlinear  effects  induced  by  the  excitation,  as  well  as 
temporal  walk-off  due  to  GVM.  We  implement  Fresnel  dif¬ 
fraction  of  the  probe  field  from  the  back  of  the  sample  to 
the  quad-cell  detector.  The  BD  signal  is  then  calculated  as 
the  difference  in  energy  detected  on  the  left-half  side  of  the 
quad-cell,  minus  that  on  the  right-half  side  [35]. 

In  typical  excite— probe  experiments  where  GVM  is 
neglected,  the  temporal  resolution  is  limited  by  the  width 
of  the  cross  correlation  of  the  excitation  and  probe  pulses. 
GVM  causes  the  two  pulses  to  overlap  at  different  distances 
within  the  sample  at  different  time  delays,  which  results  in 
a  broadening  and  smoothing  of  the  measured  signal  as  a  func¬ 
tion  of  delay  [50] .  The  temporal  resolution  of  a  given  experi¬ 
ment  is  no  longer  only  a  parameter  of  the  apparatus,  but 
depends  on  the  dispersion  of  the  material  under  investigation. 
One  consequence  of  this  broadening  is  a  reduction  of  the  tem¬ 
poral  resolution  of  the  noninstantaneous  response. 


4.  TIME  DOMAIN  EXPERIMENTS 

BD  is  used  to  measure  the  magnitude,  time  constants,  and 
polarization  dependence  of  each  nonlinear  mechanism  that 
contributes  to  the  third-order  response.  The  experimental 
setup  consists  of  slight  modifications  to  a  standard  excite- 
probe  experiment  (Fig.  J_)  and  is  the  same  as  that  reported 
in  [35] .  A  commercial  Tiisapphire  chirped-pulse-amplifier  sys¬ 
tem  (Coherent  Legend  Duo+)  with  12  mj  output  energy  at 
800  nm,  with  a  bandwidth  of  28  nm  at  a  repetition  rate  of 
1  kHz,  is  used  to  produce  the  excitation  pulse  of  duration 
50  fs  (temporal  durations  are  given  in  FWHM  unless  other¬ 
wise  specified),  as  determined  from  autocorrelation  measure¬ 
ments.  The  probe  pulse  is  generated  by  splitting  off  a 
portion  of  the  excitation  pulse  and  focusing  it  into  a  1  cm 
cuvette  filled  with  water  to  create  a  white -light  continuum, 
from  which  a  central  wavelength  of  650  nm  is  filtered  out 
using  a  narrow  bandpass  interference  filter.  Both  the  pulses 
have  Gaussian  spatial  profiles.  The  pulse  energy  of  each  beam 
is  controlled  via  a  combination  of  half-wave  plates  and  polar¬ 
izers.  The  probe  beam  is  focused  to  a  beam  waist  of 
w  =  38  pm(HWl/e2M)  at  the  sample,  and  the  excitation 
is  focused  tow,  =  170  pm.  The  excitation  beam  is  then  hori¬ 
zontally  displaced  from  the  probe  beam  by  wj 2  to  maximize 
the  deflection  angle.  The  crossing  angle  between  the  beam  is 
kept  small,  <2°,  to  prevent  spatial  walk-off  within  the 
sample.  The  deflection  of  the  probe  is  measured  by  an 
OSI  QD50-0-SD  quad-segmented  silicon  photodiode  posi¬ 
tioned  17  cm  behind  the  sample.  The  quad-cell  detector  is 
initially  centered  on  the  probe  beam.  The  excitation  induces 
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a  deflection  in  the  probe,  causing  it  to  shift  on  the  detector  and 
yield  a  difference  in  energy  on  the  left  and  right  segments, 
A E  =  -  iffight-  The  total  energy  falling  on  all  four  seg¬ 

ments,  E,  is  measured  simultaneously,  and  serves  both  for  nor¬ 
malization  and  as  a  monitor  for  nonlinear  absorption,  which 
we  do  not  observe.  Both  A E  and  E  are  lock-in  detected  at  the 
286  Hz  chopping  frequency  of  the  excitation  beam,  and  are 
measured  as  a  function  of  delay  between  the  two  pulses. 
For  a  sample  with  only  instantaneous  third-order  NLR, 
A  E/E  is  proportional  to  the  cross  correlation  of  the  pulses. 
The  pulse  width  of  the  probe  is  determined  from  a  BD  mea¬ 
surement  of  fused  silica  to  be  158  fs.  All  measurements  on  CS2 
(Sigma-Aldrich,  270660,  >99.9%)  are  conducted  using  1  mm 
path  length,  fused  silica  cuvettes.  Using  the  refractive  index 
measurements  of  [461,  we  find  p  =  4.7. 

The  refractive  index  change  seen  by  the  probe  depends  on 
the  relative  angle  6  of  its  linear  polarization  with  respect  to  that 
of  the  excitation.  The  induced  index  change  may  be 
described  by  the  index  ellipsoid  [381,  which  for  small  changes 
simplifies  to 

A n(8)  =  A«||  cos2(0)  +  Awysin2(0).  (12) 

The  relation  between  A ri\\  and  A n±  depends  on  the 
polarization  dependence  of  the  responsible  mechanisms.  For 
third-order  nonlinearities  with  isotropic  symmetry,  A = 
A«j|°/3,  and  for  reorientational,  A«J  =  -A«f|e/2.  We  may 
therefore  write  the  polarization  dependence  of  the  refractive 
index  change  as  [35] 

A n(0)  =  A«|°^cos2(0)  +  ^  sin2(6»)^ 

+  A»jf  ^cos2(8)  -  ^  sin2 (0)^ .  (13) 

The  reorientational  term  is  eliminated  at  the  so  called 
“magic  angle”  ( 6  =  arctan  \[2  «  54.7°),  leaving  only  the  iso¬ 
tropic  component  [6].  BD  measurements  are  conducted  at 
three  polarization  angles  of  the  probe  beam,  with  respect  to 
the  excitation  {8  =  0°,  parallel;  90°,  perpendicular;  and 
54.7°,  magic  angle),  which  is  adjusted  using  a  half-wave  plate 
that  is  calibrated  with  a  polarizer.  For  theses  angles 

A«(0°)  =  A«|°  +  A 
A«(90°)  =^A«f  -l~Anp 

Aw(54.7°)  =  -  Aw|°.  (14) 

The  BD  data  is  fit  using  Eq.  (10).  We  begin  at  the  magic 
angle,  since  it  isolates  the  isotropic  components  and  allows 
unambiguous  determination  of  their  magnitudes  and  time 
constants.  As  shown  in  Fig.  2,  the  signal  at  the  magic 
angle  is  dominated  by  the  bound-electronic  mechanism  but 
also  exhibits  a  small  noninstantaneous  contribution  consistent 
with  other  observations  [6,7,  j_0].  We  attribute  this 


Fig.  2.  Measured  (circles)  and  fit  (curves)  BD  signals  for  parallel 
(black),  magic  angle  (blue),  and  perpendicular  (red)  polarizations;  inset 
shows  the  resultant  noninstantaneous  responses  of  collision  (blue),  dif¬ 
fusive  reorientation  (red),  libration  (green),  and  their  sum  (black). 


noninstantaneous  contribution  to  the  collision-induced 
nuclear  mechanism  [6],  due  to  its  symmetry  and  fit  with 
Eq.  (8).  The  magnitudes  are  scaled  by  9/5  to  give  A»jj® 
[Eq.  (14)],  which  is  then  applied  to  the  measurements  with 
parallel  and  perpendicular  polarizations,  which  show  two 
additional  noninstantaneous  components  corresponding  to 
the  libration  and  diffusive  reorientation  mechanisms.  These 
are  fit  with  Eqs.  (5)  and  (6),  respectively.  We  do  not  resolve 
the  rise  times  for  the  collision-induced  and  diffusive  reorien¬ 
tation  mechanisms,  and  therefore  follow  [1],  and  assume  them 
to  be  equal  since  inertial  delay  should  influence  each  similarly. 
The  best  fit  of  both  parallel  and  perpendicular  polarization 
data  sets  is  used  to  determine  the  final  values  given  in  Table  L 
Errors  in  n2i m  are  estimated  from  the  uncertainty  in  irradiance 
(~20%)  and  the  signal-to-noise  ratio.  Those  for  the  remaining 
parameters  are  determined  in  the  fitting  procedure  by  varying  a 
single  component  until  the  curves  no  longer  match  the  data. 
Our  results  agree  well  with  previous  measurements  using 
similar  models  but  here  include  absolute  magnitudes  as 
well  as  response  times  [1,4,5,9,34].  Using  Eq.  (2),  x'Xxxx  = 
(1.8±0.4)x  10-21  m2/V2  [(1.3±0.3)x  10-13  esu],  and 

X%]xy  =Xxyyx  =  Xmcyy  =  CtEL/3)  =  (0.6  ±  0. 1)  X  10-21  m2/V2 
[(0.4±0.1)x  10-13  esu]. 

Using  these  results  we  are  then  able  to  predict  the  results  of 
a  DFWM  experiment.  This  was  conducted  using  another 
Tirsapphire  chirped-pulse-amplifier  (Clark-MXR  CPA  2010) 
with  1  mj  output  energy  at  a  repetition  rate  of  1  kHz,  which 
pumps  an  optical  parametric  amplifier  (Light  Conversion, 
TOPAS-C)  tuned  to  700  nm.  A  prism  compressor  (SF10) 


Table  1.  Fit  Parameters  of  Third-Order  Response  of  CS2a 


Mechanism 

n2,m 

Cv„  (fs) 

T/,m  (fs) 

Symmetry 

Electronic 

2.0  ±  0.4 

Instantaneous 

iso 

Collision 

1.0  ±  0.2 

150  ±  50 

140  ±  50 

iso 

Libration 

7.6  ±  1.5 

b 

450  ±  100 

re 

Diffusive 

18  ±  3 

150  ±  50 

1610  ±  50 

re 

a^2,m  are  given  in  units  of  10  19  m2/W. 
bco0  =  8.5  ±  1.0  ps-1  and  a  =  5  ±  1  ps-1. 
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Fig.  3.  Normalized  DFWM  signal  of  CS2  (circles),  and  calculation 
(curve)  using  the  response  function  model  values  of  Table  1_;  inset  shows 
logarithmic  scale. 


compresses  the  pulse,  which  is  then  split  by  a  50/50  beam¬ 
splitter,  and  the  reflected  pulse  is  sent  to  a  delay  line.  The 
two  beams  are  overlapped  onto  a  transmissive  diffraction 
grating.  The  ±1  diffracted  orders  of  the  top  beam  act  as 
the  excitation,  and  the  +1  order  of  the  bottom  beam  acts 
as  a  probe.  The  two  excitation  beams  are  focused  and  over¬ 
lapped  in  the  CS2  and  generate  a  refractive  index  grating 
off  of  which  the  probe  beam  diffracts  [9,J_0].  The  measured 
signal  is  proportional  to  the  diffraction  efficiency,  which  is 
related  to  the  square  of  the  induced  refractive  index  change. 
The  prism  compressor  is  adjusted  to  yield  the  narrowest  signal 
in  delay.  Figure  3  shows  the  normalized  DFWM  signal  versus 
delay  along  with  a  simulation  of  the  experiment  using  the  re¬ 
sponse  function  with  the  values  in  Table  J_  for  88  and  42  fs 
pulse  widths  of  the  excitation  and  probe  pulses,  respectively. 
The  calculation  shows  excellent  agreement  with  the  experi¬ 
mental  data,  validating  the  results  obtained  with  the  BD  tech¬ 
nique,  and  demonstrates  our  ability  to  apply  the  measured 
response  function  to  predict  the  outcome  of  other  NLO 
experiments  on  CS2. 


filtered  to  produce  a  Gaussian  beam.  A  small  portion 
(~10%)  of  the  beam  is  split  off  to  a  reference  detector  to  win¬ 
dow  the  pulse  energy,  which  is  controlled  by  a  half-wave  plate 
and  polarizer.  A  quarter-wave  plate  is  inserted  before  the  focus¬ 
ing  lens  for  circular  polarization  measurements.  A  beam  splitter 
after  the  sample  allows  both  open  aperture  (OA)  and  closed 
aperture  (CA)  Z-scans,  which  are  sensitive  to  NLA  and 
NLR,  respectively,  to  be  performed  simultaneously.  A 
mode-locked  Nd:YAG  laser  system  (EKSPLA  PL2143)  oper¬ 
ating  at  10  Fiz  is  also  used  to  measure  n2,eff  with  ps  pulses.  In 
addition  to  the  1064  nm  laser  output  itself,  which  had  a  pulse 
width  of  25  ±  1  ps,  it  is  also  frequency  tripled  and  used  to 
pump  an  OPG/A,  which  is  tuned  to  700  nm  with  a  pulse 
width  of  8  ±  1  ps.  For  all  measurements  the  thin  sample 
approximation  is  valid  [16,511- 

OA  Z-scan  measurements  of  2PA  in  ZnSe  for  700  nm  (with 
2PA  coefficient  a2  =  5.6  x  10~u  m/W  [52]),  and  CdSe  for 
1064  nm  (a2  =  23  x  10-11  m/W  [52]),  and  CA  Z-scan  mea¬ 
surements  of  fused  silica  (n2,e/  =  0.25  x  10~19  m2/W  [53]) 
are  used  to  determine  the  spot  size  and  pulse  width  by  fitting 
with  the  standard  Z-scan  analysis  [16] .  Determining  the  pulse 
width  and  spot  size  in  this  manner  minimizes  the  relative  error 
between  Z-scans  at  different  pulse  widths.  Z-scans  of  CS2  are 
performed  at  multiple  energies  to  ensure  the  linear  dependence 
of  the  peak-to-valley  transmission  change  with  irradiance,  and 
the  best  fit  to  all  curves  is  used  to  determine  the  value  of  n2?es 
(see  Supplement  1). 

The  Z-scan  technique  is  sensitive  to  changes  of  the  materi¬ 
al’s  refractive  index,  which  are  only  proportional  to  «2  (,/,  in  the 
case  of  an  instantaneous  response.  Therefore,  an  irradiance- 
weighted  time  average  has  to  be  performed  in  order  to  map 
the  BD  results  on  the  pulse  width  dependent  n2eg  values  de¬ 
termined  from  the  Z-scan  measurements.  The  time  averaged 
index  change  over  the  duration  of  the  pulse  is  [16] 


(A«(r)} 


f  I(t)dt 


(15) 


For  a  response  assumed  to  be  instantaneous 


5.  Z-SCAN  EXPERIMENTS 

CS2  is  often  used  as  a  reference  material  for  Z-scan  experiments 
despite  its  pulse  width  dependent  NLR.  From  the  measured 
third-order  response  function  the  effective  nonlinear  refractive 
index,  defined  by  A n  =  n2egl,  can  be  calculated  and  com¬ 
pared  to  Z-scan  experiments.  Defining  n2eg  in  this  way  is 
commonly  done  but  does  not  address  the  dynamics  of  the 
physical  processes  that  cause  the  index  change  and 
results  in  pulse  width  dependent  values.  To  check  our  predic¬ 
tions  of  »2jeff,  we  conduct  Z-scans  at  both  700  and  1064  nm, 
over  a  broad  range  of  pulse  widths  by  using  a  prism  compres¬ 
sor/stretcher.  The  same  commercial  Ti:sapphire  laser  system 
used  for  the  DFWM  experiments  is  also  used  here.  The  exper¬ 
imental  setup  has  been  described  previously  [16].  The  prism 
compressor  is  used  to  vary  the  pulse  duration  from  32  fs 
(corresponding  to  the  bandwidth  limit)  to  2.3  ps  (limited 
by  available  space  in  the  setup).  The  beam  is  then  spatially 


(A  n(t))  =  nXcS 


f  1 2  (t)dt 
f  I{t)dt 


(16) 


Equating  Eqs.  (15)  and  (16),  using  Eq.  (J_)  for  A n(t),  and 
solving  for  n2  tg,  yields 


^2, eft  n'l,d  + 


S I {t)  f  R(t-t')I(t')dt'dt 
f  I2(t)dt 


(17) 


Figure  4  shows  the  measured  values  of  n2  eff  from  Z-scans 
using  linear  polarization  to  those  calculated  using  Eq.  (17). 
Measurements  at  both  700  and  1064  nm  overlap,  demonstrat¬ 
ing  the  negligible  dispersion  of  the  noninstantaneous 
contribution  to  the  NLR.  In  the  limit  of  short  pulses, 
n2,eff  =  n2 ,eb  since  only  the  bound-electronic  response  con¬ 
tributes.  As  the  pulse  width  increases,  the  noninstantaneous 
contributions  add  to  n2  tg,  which  plateaus  for  pulses  longer 
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Fig.  4.  Comparison  ofZ-scan  measurements  (see  Supplement  1)  using 
the  Ti:sapphire  (dosed)  and  Nd:YAG  laser  system  (open),  at  both 
700  nm  (black)  and  1064  nm  (green),  and  calculation  using 
Eq.  (17),  red  curve,  of  of  CS2  versus  pulse  width;  shaded  region 
represents  errors  in  the  response  function  from  Table  1. 


than  ~10  ps.  In  this  long-pulse  limit  n2,eg  *s  given  by  the  sum 
of  each  component’s  magnitude: 

^2, eft t long  ^  T2, m-  0 

m 


Our  results  agree  with  the  BD  predictions,  using  the  values 
of  Table  J_,  and  with  other  Z-scan  measurements  of  »2,eff 
[16-21,54].  Z-scan  measurements  are  also  performed  using 
circular  polarization.  The  ratio  of  n2?eg  f°r  linear  and  circular 
polarizations,  W2eff/W2eff>  *s  an  indication  of  the  polarization 
dependence  of  the  underlying  mechanisms.  For  materials  pos¬ 
sessing  isotropic  symmetry,  this  ratio  is  1.5,  while  for  reorienta- 
tional  mechanisms  it  is  4  [381.  In  CS2,  »2"ff /«“reCfj-  depends  on 
the  pulse  duration.  By  using  the  symmetry  properties  of  the 
response  function,  we  find  A»circ  =  A»y°/1.5  +  A«j|e/4. 
Figure  5  shows  excellent  agreement  between  the  measured 
values  «2neff /  w2lreCff  versus  pulse  width  and  that  predicted  by 


Fig.  5.  Comparison  of  «2dt/”zSr  versus  pulse  width  between  Z-scan 
measurements  with  both  Ti:sappire  (closed  circles)  and  Nd:YAG  (open 
circle)  laser  systems  at  700  nm  and  calculated  (red  curve);  shaded  region 
represents  only  relative  errors  that  contribute  to  uncertainty. 


the  third-order  response  function  in  Table  L  In  the  short  pulse 
limit,  «2eff/ w2eff  =1-5  because  it  is  dominated  by  the  iso¬ 
tropic  bound-electronic  response,  while  in  the  long  pulse  limit 
this  ratio  becomes 


Jin 

n2,eff 


'2,eff 


Mlin 


'2,effllong 


long 


(»2 ,'l+nZc)  I  + 

“V  4 


1.5 


(19) 


where  the  numerator  is  given  by  Eq.  (18),  and  the  denomina¬ 
tor  shows  the  weighting  of  each  contribution  according  to  its 
polarization  dependence.  Using  the  results  in  Table  J_,  for  long 
pulse  widths  «2”ff/«2eff  =  3-4. 

To  study  the  dispersion  of  the  bound-electronic  component 
of  the  third-order  response,  we  conduct  Z-scans  at  several 
wavelengths  from  390  to  1550  nm,  at  the  shortest  pulse  widths 
achievable.  Depending  on  the  wavelength,  the  minimum  pulse 
width  varies  from  32  to  165  fs.  Over  this  range,  the  nonlin¬ 
earity  of  CS2  is  dominated  by  the  bound-electronic  response, 
but  the  noninstantaneous  component  can  still  contribute  sig¬ 
nificantly,  and  can  increase  n2:eg  up  to  ~2x.  However,  if  we 
assume  the  nuclear  response  is  independent  of  wavelength,  we 
can  calculate  its  contribution  for  each  particular  pulse  width, 
and  simply  subtract  it  out.  That  is,  since  we  know  R(t)  from 
the  BD  measurements,  and  we  measure  «2>e g-  with  a  Z-scan,  we 
can  calculate  n2ei  using  Eq.  (17).  To  verify  our  assumption,  we 
perform  Z-scans  using  long  pulses,  both  with  the  ps  system  as 
well  as  a  10  Hz  Q-switched  Nd:YAG  (Continuum  PL9010)  ns 
system,  pumping  an  optical  parametric  oscillator  (Continuum 
Sunlite  EX),  that  produces  pulses  from  2.5  to  16.5  ns,  depend¬ 
ing  on  the  wavelength.  For  ns  pulses,  electrostriction  can  also 
contribute  to  the  closed  aperture  Z-scan  signal,  and  analysis  is 
done  to  extract  its  contribution  [54],  Figure  6(a)  shows 
»2,eff  [long  *s  independent  of  wavelength  to  within  the  measure¬ 
ment  uncertainty.  Since  the  nuclear  components  dominate  for 
long  pulses,  we  conclude  that  they  have  negligible  dispersion. 
Therefore,  for  the  short  pulse  measurements  we  may  subtract 
the  nuclear  contribution  to  n2,es  from  the  value  measured  via 
Z-scan,  and  isolate  only  the  bound-electronic  contribution. 

Figure  6(b)  shows  both  the  dispersion  of  n2>ei  with  the  nu¬ 
clear  contribution  removed  and  the  spectrum  of  a2.  A  quan¬ 
tum  mechanical  perturbation  approach,  known  as  the  “sum 
over  states”  (SOS)  model,  which  considers  all  electronic  dipole 
allowed  transitions  between  various  states,  is  used  to  determine 
the  NLO  response  [55].  In  this  way  n2:ei  and  a2  can  be  calcu¬ 
lated  from  knowledge  of  the  energies  of  the  electronic  states 
and  their  respective  transition  dipole  moments  ji.  Here,  we 
use  only  three  states  in  the  model  [56,57]:  the  ground  state, 
a  one-photon  allowed  state  e,  and  a  2PA  final  state  e' .  Due  to 
the  symmetry  of  the  molecule,  the  states  g  and  e  have  the  same 
parity  and  thus  [ige’  =  0.  The  results  are  shown  in  Table  2. 
Parameters  for  state  e  were  found  by  measuring  the  linear  ab¬ 
sorption  spectrum  of  a  0.445  mM  solution  of  CS2  in  meth¬ 
anol.  The  reported  transition  dipole  moments  are  vacuum 
values,  i.e.,  they  have  been  corrected  for  local  field  effects 
[38] .  The  corresponding  spectral  shape  of  n2>ei  agrees  well  with 
the  experimental  data  both  in  peak  position  and  trend  of 
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Fig.  6.  (a)  Z-scan  measurements  of  dispersion  of  it2,efifli0ng  f°r  (°Pen 

circles)  13  to  20  ps  and  (closed  squares)  2.5  to  17  ns,  with  the  electro- 
striction  contribution  subtracted;  (b)  Z-scan  measurements  of  NLR 
(black  circles)  for  fs  pulses  with  noninstantaneous  component  subtracted 
and  a2  (blue  triangles);  curves  represent  the  SOS  model  fit  for  2PA  (blue) 
and  n2:ei  (black),  which  has  been  multiplied  by  a  factor  of  2. 


Table  2. 
CS2 

Fit  Parameters  for  SOS  Model  of  n2,ei  and  a2  of 

State 

Energy  (eV) 

HWHM  (eV) 

A  (D) 

e 

t 

e 

6.00  ±  0.01 
5.95  ±  0.05 

0.17  ±  0.01 

0.40  ±  0.05 

Hge  =  3.6  ±  0.2 
A„'  =  8.7  ±  0.9 

variation.  However,  it  underestimates  the  ma 
tor  of  2,  which  may  be  a  consequence  of  the 

gnitude  by  a  fac- 
simplicity  of  the 

model.  Using  these  parameters,  the  SOS  model  may  be  applied 
to  predict  the  values  of/*3*  for  other  optical  frequencies  (e.g., 
third  harmonic  generation  (THG)  [15])- 

6.  CONCLUSION 

By  the  application  of  several  measurement  techniques,  we  have 
experimentally  investigated  the  third-order  NLO  response  of 
CS2  to  give  a  complete  picture  of  its  temporal,  polarization, 
and  spectral  dependence.  Using  the  recently  developed  BD 
technique  [35],  we  have  measured  the  absolute  magnitude 
of  the  third-order  response  function.  We  have  expanded  the 
analysis  of  the  BD  technique  to  include  a  noninstantaneous 
component  of  NLR  and  GVM.  A  model  of  instantaneous 
bound-electronic  and  three  non-instantaneous  nuclear  re¬ 
sponses,  i.e.,  collision,  libration,  and  diffusive  reorientation, 
is  found  to  give  excellent  agreement  with  the  BD  experimental 
data.  The  polarization  dependence  of  each  mechanism  is  used 


to  determine  their  respective  tensor  symmetries.  Calculations 
of  the  results  of  DFWM  using  the  BD  generated  response 
function  give  excellent  agreement  with  experimental  results. 

Based  on  the  measured  response  function,  we  are  able  to 
predict  the  pulse  width  dependent  n2>eg ,  which  agrees  very  well 
with  Z-scan  measurements  over  a  range  from  32  fs  to  17  ns. 
For  long  pulses  slower  effects  also  contribute  to  the  refractive 
index  change,  including  electrostriction  (ns)  and  thermal  ef¬ 
fects  (ms),  which  have  been  well-characterized  previously 
[54,58].  With  the  polarization  dependence  of  the  various  re¬ 
sponse  mechanisms,  predicted  results  of  Z-scans,  conducted 
with  both  linear  and  circular  polarizations,  also  give  excellent 
agreement  with  experimental  results.  Since  the  dispersion  of 
the  nuclear  mechanisms  is  shown  to  be  negligible,  the 
dispersion  of  the  bound-electronic  response  is  isolated  in 
short-pulse  Z-scan  measurements.  The  dispersion  of  n2e/ 
and  the  spectrum  of  2PA  across  a  broad  wavelength  range 
(from  390  to  1550  nm)  is  presented.  This  allows  analysis 
in  terms  of  a  simplified  quantum  mechanical  model  that 
matches  the  spectral  shape. 

We  are  therefore  able  to  predict  the  nonlinearity  of  CS2 
over  a  broad  range  of  experimental  parameters.  This  paper  will 
prove  useful  for  those  making  use  of  CS2  for  a  multitude  of 
NLO  applications  including  as  a  reference  material  for  NLO 
spectroscopy.  The  experimental  approach  demonstrated  here 
may  be  applied  to  a  wide  variety  of  other  materials  for  robust 
quantitative  determination  of  the  temporal,  spectral,  and 
polarization  dependence  of  their  NLO  response. 
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